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Many of the exotic properties proposed to occur in graphene rely on the possibility of increasing 
the spin orbit coupling (SOC). By combining analytical and numerical tight binding calculations, 
in this work we study the SOC induced by heavy adatoms with active electrons living in p orbitals. 
Depending on the position of the adatoms on graphene different kinds of SOC appear. Adatoms 
located in hollow position induce spin conserving intrinsic like SOC whereas a random distribution 
of adatoms induces a spin flipping Rashba like SOC. The induced SOC is linearly proportional to 
the adatoms concentration, indicating the inexistent interference effects between different adatoms. 
By computing the Hall conductivity we have proved the stability of the topological quantum Hall 
phases created by the adatoms against inhomogeneous spin orbit coupling . For the case of Pb 
adatoms, we find that a concentration of 0.1 adatom per carbon atom generates SOC’s of the order 
of ~40mey. 


I. INTRODUCTION 

The research on graphene, a two-dimensional crystal 
of carbon atoms, has driven to the discover of a large 
number of interesting electrical, magnetic, mechanical 
and optical properties[Tl[2]. The small atomic number of 
carbon implies that electrical carriers in graphene have 
a extremely weak spin-orbit (SO) coupling[3l [4]. This 
property in combination with the large graphene elec¬ 
tron mobility makes graphene a very good candidate for 
using in spintronics[5l E]. 

On the other hand some proposasl of exotic topologi¬ 
cal phases in graphene rely on the possibility of increas¬ 
ing the SOC. Because of the graphene lattice symmetry, 
there are two types of spin orbit coupling in graphene, 
intrinsic-\\k.e^ where the ^-component of the electron spin 
is a good quantum number and RashhaASke which mixes 
spins and appears in absence of mirror symmetry [ 7 ]. In 
graphene, the intrinsic SOC opens a gap and the sys¬ 
tem becomes a quantum spin Hall insulator, with gap¬ 
less edge states able to transport spin and charge [H |9]. 
Non-trivial topological phases may also occur in multi¬ 
layer graphene m- Quantum anomalous Hall effect was 
predicted to occur in bilayer and monolayer graphene in 
presence of Rashba SOC and an exchange field or mag¬ 
netic impurities pTl [T2]. The experimental realization 
of these topological phases requires a large SOC, and 
therefore there is a big interest [T3HT9] in increasing the 
SOC and clear the way to the study of exotic topological 
phases in graphene. Experimental reports on enhance¬ 
ment of SOC in graphene by weak hydrogenation [17], 
gold hybridization [20] or proximity with WS 2 m, indi¬ 
cate that it is possible to increase the SOC in almost three 
orders of magnitude. Recently, it has been reported that 
graphene grown on Cu shows a SO splitting around 20 
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meV|21|. Intercalation of An atoms in graphene grown 
on Ni produces a SO splitting of near 100meV[20]. Sim¬ 
ilarly, the intercalation of Pb atoms in graphene grown 
on a iridium substrate seems to produce a giant SOCfTSj. 
Theoretically, it has been proposed that heavy adatoms 
with partially filled p-shells, deposited on symmetric po¬ 
sitions of the graphene lattice, could induce large intrinsic 
SOC[22]. 



FIG. 1. (Color online)Schematic representation of the effec¬ 
tive hopping between carbon atoms induced by an adatom 
with active electrons in outer p-shells. Vertical, H and T, 
lines indicate the hollow and top positions respectively 

In this work we study the SOC induced by heavy 
adatoms with active electrons living in p-orbitals, in par¬ 
ticular we consider Pb atoms, that have been proved to 
induce large SO effects in grapheneflS]. 

The physical picture is the following, the tunneling 
of an electron between two carbon atoms through the 
adatoms p-orbitals opens new channels for hopping in 
graphene. The SOC between the adatom p orbitals 
makes that the new tunneling channels can conserve the 
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spin inducing a intrinsic SOC or can flip the spin inducing 
a Rashba like SOC. 

By combining analytical calculations, perturbation 
theory and tight-binding based numerical simulations, we 
study the type of SO coupling induced by adatoms resid¬ 
ing in different positions of the graphene unit cell. In 
addition, we study how a finite density of adatoms, in 
different distributions, affects the induced SO couplings. 
The main conclusions of our work are the following, 

i) adatoms located in hollow positions, see Fig{^ induce 
intrinsic SOC. A finite density of adatoms in hollow po¬ 
sitions opens an energy gap at the Dirac points, that 
increases linearly with the adatom concentration. This 
gapped phase is a quantum spin Hall state. The simula¬ 
tions indicate that, even for high adatom coverage, there 
are not interference effects between the adatoms and the 
gap only depends on adatom density. In the case of Pb 
atoms we find gaps of the order of 50meV for a concen¬ 
tration of 0.1 adatom per carbon. 

ii) for adatoms placed in top positions, see Fig{^ the tun¬ 
neling from graphene to the adatom and back, induces a 
Rashba like spin flip hopping between the underneath C 
atom and its first neighbors and an intrinsic like spin con¬ 
serving second neighbors tunneling between the carbons 
surrounding the underneath carbon atom. The intrinsic 
like SOC induced by adatoms in top positions has oppo¬ 
site sign than the induced by adatoms in hollow geometry. 
The Rashba SOC has the same sign independently of the 
sublattice of the underneath carbon. 

iii) A finite density of adatoms randomly distributed on 
graphene, induces a finite Rashba SOC linearly propor¬ 
tional to the density of adatoms. For a random distri¬ 
bution of adatoms, the resulting intrinsic like SOC van¬ 
ishes, because contributions from different locations of 
the adatoms have opposite signs. Similar results are ob¬ 
tained when the adatoms form an array commensurate 
with a large graphene supercell. By computing the Hall 
conductivity, we have obtained that a random distribu¬ 
tion of adatoms on graphene in presence of an exchange 
field, is an anomalous quantum Hall system. When the 
adatom is Pb, we obtain that the Rashba SO coupling 
can be as larger as 35meV for a concentration of 0.1 Pb 
per carbon atom. 

The rest of the paper is organized in the following 
way, in Section H we introduce the graphene and adatom 
Hamiltonians, and in Section HI we describe the hop¬ 
ping between graphene carbon atoms and the adatom 
p-orbitals. In Section IV we present the perturbation 
theory for describing the adatom mediated effective hop¬ 
ping between carbon atoms. The knowledge of the effec¬ 
tive hopping between carbon atoms allow us to obtain, 
in Section V, analytical expressions for SOC induced by 
adatoms located in top and hollow positions. Section VI 
turns to present tight-binding based numerical simula¬ 
tions for studying the effect that a random or commen¬ 
surate distribution of adatoms have on the induced SOC. 
In Section VH we calculate the topological properties of 
graphene doped with adatoms. We close the paper with 


a summary of the results. 

II. PRELIMINARIES. 

A. Graphene Hamiltonian 

In graphene, carbon atoms crystallize in a triangu¬ 
lar lattice of primitive translation vectors a=(0,a) and 
b=(v^/2, l/2)a, where a=2.46A is the lattice constant. 
The positions of the triangular lattices are R^. There are 
two atoms per unit cell located at positions dA=(0,0) 
and dB = (a/v/3, 0), that define sublattices A and B in 
graphene. Covalent sp 2 bonds between carbon atoms 
stabilize this honeycomb lattice whereas the tunneling 
between pz orbitals is the origin of the low energy active 
conduction and valence 7r-bands. The band structure is 
rather well described by a tight-binding model with hop¬ 
ping t 2.7eV between first neighbors carbon Pz orbitals, 

Ho = -t y] {\Z,iA,cr >< Z,iB,cr\ + H.c.) (1) 

<iA,iB>,cr 

Here the sum runs over first neighbors pairs and 
\Z,ia^cF > represents the wavefunction of an electron at 
position Ri+do, occupying a carbon pz orbital with z- 
component of the spin a. The energy of the Pz carbon 
orbital is chosen as the zero of energies. 

In graphene the intrinsic SOC has a chiral structure of 
the form, 

<ia,ja>,cr 

( 2 ) 

where the sum runs over second nearest neighbors carbon 
atoms and is a unit vector parallel to Ri+dc,. Note 
that the intrinsic SOC conserves spin and it is not asso¬ 
ciated with broken mirror symmetry. On the contrary, 
Rashba SOC appears because of broken mirror symme¬ 
try, in particular due to the substrate, and induces a 
coupling between first neighbors with opposite spin of 
the form, 

Hso=i^^ V ii'^'^'^iAjB)z\Z,iA,':r><Z,jB,<y'\ +H.c.) 

( 3 ) 

here Ui^j^ is a unit vector parallel to Rj + d/3 — R^ — 
and a the electron spin Pauli matrices. 

In absence of SO couplings the conduction and valence 
bands touch at two inequivalent points of the Brillouin 
zone K=(0,1^) and K'=(0, —|^) which are the cele¬ 
brated Dirac points. Near these points the low energy 
physics is described by the Dirac equation 

Ho = hvF{kxO-o O + skyao (g) Ty) (4) 

here the moment k is measured with respect the Dirac 
points, 5=1 and -1 stands for K and K' respectively, and 
T are the Pauli matrices acting on the spinor defined 
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by the amplitude of the wave function on sublattices A 
and B. In the previous equation cfq represents the unity 
matrix in the spin sector. The Fermi velocity is related 
with the hopping trough the relation hvF=^ta. In this 
continuum approximation the SO terms get the form, 

HgQ = (Tz (8) Tz (5) 

3 

^SO = crx®Ty-ay® Tx) . ( 6 ) 

B. Adatom Hamiltonian. 

We consider heavy atoms with the active electrons liv¬ 
ing in p orbitals. The Hamiltonian describing the elec¬ 
trons in the adatom contains a spin-orbit coupling part 
and a crystal field Hcf term. In the basis{|pa, \py t> 
5 \Pz t >5 \Px \Py \Pz i>} the Hamiltonian reads, 


Hp = 

A^oL • 

a + Hcf - 
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Here Aso is the spin-orbit coupling parameter, L and 
cr are the usual angular momentum and spin operators. 
Non spherical effects occurring in the geometry produce 
a crystal field that splits the energies of the p orbitals. 
For adatoms deposited on planar graphene, we expect 
that ex=ey ^ tz- 



FIG. 2. (Color online)Geometry of and adatom located at a 
distance h of the graphene layer. The angles 0 and 0 and the 
distance d define the spherical coordinates of the adatom with 
respect to a carbon atom. 


III. TUNNELING BETWEEN AN ADATOM 
AND A CARBON Z ORBITAL. 

We consider an adatom placed at position r=(x,^, h), 
where h is the vertical distance between graphene and 
the adatoms, Fig§ The tunneling amplitudes between 
an carbon orbital located at position Ri = (Xi, T^, 0 ), and 
the adatom py and Pz orbitals are, 

Z^ (j^T^px-, O' ^ — — cos 0 sin 2 ^(I^pcr(d) I^p 7 r(d)) 

Z, i, (T|T|py, (T ^ — — sin 0 sin 2 d(I^pcr(d) I^p 7 r(d)) 

< Z,i,o\f\pz,o > = cos^ 0Vpp^{d) 0Vppa{d) ( 8 ) 

where d=tan“^ , ^ , d)=tan“^ , and 

Vppa- and VppT^ are the Slater-Koster hopping parameters 
between the Pz graphene orbital and the heavy adatom p- 
orbitals. The hopping parameters depend on the distance 
d= between the atoms, that 

we parametrize in the form, Ippcr( 7 r) (d)=Ippcr( 7 r) (d = 
with /3=3[23] and Ippcr( 7 r)(d) obtained from 
density functional calculations[THj. Note that in the tun¬ 
neling process the carrier spin is conserved. Because 
of the symmetry of the p-orbitals, the hopping ampli¬ 
tude between a carbon p^-orbital in graphene and the 
Px and Py adatom orbitals, have opposite sign depending 
weather the adatom is deposited on top or bottom of the 
graphene sheet. On the contrary, the hopping between 
z-orbitals is independent of the position of the adatom 
with respect the graphene layer. 

IV. TUNNELING BETWEEN CARBON ATOMS 
MEDIATED BY AN ADATOM. 

The spin-orbit coupling in the adatom located at 
r=(x,^,d), allows an extra path for tunneling between 
two carbon atoms located at R^ and Hj with spin o and 
o' respectively. In second order perturbation theory a 
single adatom produces a coupling between the carbon 
atoms of the form, 

s:^<Z,i,a\f\pi><pi\f\Z,j,a' > 

^ ~ 

I ^ 

here \pi > and q are the eigenfunctions and eigenvalues 
of Hamiltonian Eqj^ Because of the form of spin-orbit 
coupling in the adatom outer p shell, the induced hop¬ 
pings satisfy the following relations for i ^ j 

_ 5l< 

—CTjjfcr 

_ 

T*cr,jicr —cr,jf —cr 

7ia,ja' = lja',ia (^0) 

and for i = j 

"Jiajia Ti —cr,i —cr 

Ticrji —cr — 0 


( 11 ) 

( 12 ) 
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The adatom gives rise to two kind of SO assisted tun¬ 
neling, 

i) spin conserved tunneling events of the form 
\Z,i,a>^ \Px{y)cr>^ \Py{x)cr>^ \z,j,a> (13) 

that are pure imaginary and change sign when reversing 
spin. Following the standard notation we call it intrinsic 
spin-orbit coupling. This tunneling amplitude behaves as 
sin^(26>), it is zero when the adatom is in the graphene 
sheet, being independent on the top or bottom position 
of the adatom with respect the graphene layer. 

ii) non-conserving spin processes of the form 

\Z,i,(j>^ \p^(j>^ \Px(y) -<7>^ \Z,j,-a> 

\Z,i,a>^ \Px{y)<y>-^ \Pz -cr>^ |Z,j,-(7>(14) 

These terms get origin on the lack of mirror symmetry 
in the hopping between z-orbitals and we refer to this 
tunneling contribution as Rashba SOC. This tunneling 
amplitude behaves as sin(2^) and changes sign when the 
adatom is located on top or on bottom of the graphene 
layer. 

V. EFFECTIVE HAMILTONIANS FOR 
HOLLOW AND TOP POSITIONS 

When the adatoms are located at high symmetry 
points of the graphene lattice, it is possible to write down 
analytic expressions for the effect that the SO induces 
on the graphene low energy band structure. The pro¬ 
cedure consist in projecting the perturbation created by 
the adatom on the atomic Bloch states at the K and K' 
Dirac points, 

^ i 

= (15) 

here N is the number of of unit cells in the crystal. These 
Bloch states are the eigenstates of the Dirac Hamilto¬ 
nian, Eqj^for k=0. We assume that adatoms do not in¬ 
duce coupling between states coming from different Dirac 
cones. We have checked numerically this assumption, 
provided the adatoms do not form a periodic array with 
a reciprocal lattice vector equal to K — K'. 

A. Adatom in Hollow Position 

In the hollow geometry the adatom is located on top 
of the center of an hexagon of the graphene lattice at a 
height h, see Figl3 We consider SO induced tunneling 
up to third neighbors, tunneling between more distant 
atoms can be neglected because of the exponential de- 
creasement of the tunneling amplitude with the distance. 


Coupling between Bloch wavefunctions of different sub¬ 
lattices involves first and third neighbors hopping and 
gets the form 

< >= ^ E 

iA,jB 

(16) 

here iA {js) runs over the vertices, of sublattice A (H), of 
the hexagon surrounding the adatom. Vh represents the 
perturbation created by the adatom in hollow position. 

The coupling between Bloch states of the same sub¬ 
lattice involves second neighbors tunneling and gets the 
form, 

iAT^jA 

(17) 

Similar expression applies for < B,s,(T\^H\'^\B,s,cr' >• 

The adatom also induces diagonal selfenergies that for 
the adatom in the hollow position are equal for both 
Dirac points, spin orientation and graphene sublattices. 

In the hollow geometry and for cr'=—cr, it is possible 
to sum the six terms, Eqj^ that contribute to spin flip 
effective tunneling, and we get 

= -atR ( 18 ) 

being Ir a constant that depends on the carbon to 
adatom tunneling parameters, Eq§ Using this expres¬ 
sion and relations Eq{^ we obtain that in the hollow 
position an adatom with outer shell p orbitals does not 
induce non-conserving spin tunneling and therefore does 
not induce Rashba like SOC in graphene. 

Eor spin conserving SO induced tunneling, the sum of 
the two processes described in Eq(T^ gives a hopping, 

liaja = iCFtso siu {(j)i - (j)j) (19) 

where tgo is a constant that depends on the distance be¬ 
tween the adatom and the graphene sheet. When intro¬ 
ducing this hopping and applying the symmetries Eqj^ 
we obtain that the coupling between Bloch states of dif¬ 
ferent sublattices and same spin cancels identically. On 
the contrary the spin conserving coupling between same 
sublattice Bloch functions gets a finite value that changes 
sign when changing spin, sublattice or Dirac cone, 

< '^r,s,a\yn\'^T,s,a >= ^^VStgoCTST . ( 20 ) 

This term has the same form than the Hamiltonian Eqj^ 
and we conclude, in agreement with reference [H], that 
an heavy adatom, with electrical active p-orbital, in a 
hollow position on top of graphene induces an intrinsic 
like SOC. 

B. Adatom in Top Position 

In this geometry the adatom is located vertically on top 
of a carbon atom at a height h. This configuration priv¬ 
ileges the sublattice A of the underneath carbon atom. 
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In the top arrangement the carbon orbital is orthog¬ 
onal to the Px and Py orbitals of the adatom located on 
top of it. Therefore the processes contributing to first 
neighbors spin conserving tunneling are zero by symme- 
try, Eq|T3l However, an adatom on top of a carbon of 
a given sublattice, induces spin conserving tunneling be¬ 
tween carbon atoms of the opposite sublattice, 

3v^ 

<^S(A),S,cr|^TA(B) I^B(A),S,cr >= tso ^ (21) 

here Vta^b) represents the perturbation created by the 
adatom on top of atoms belonging to sublattice A{B). 
Therefore adatoms in top positions induce intrinsic-like 
SO coupling, although it is important to note that the 
sign of this conserving tunneling is opposite to the in¬ 
duced by an adatom in hollow position, Eqj^ 

Because of the symmetry of the orbitals, only one of 
the mechanisms described in Eq{^ contributes to spin 
flip tunneling between first neighbors, 

, ( 22 ) 

where Ir is a constant that depends on carbon to adatom 
tunneling parameters. Adding the contributions from the 
three first neighbors of the underneath C atom, we get 
the following contribution to the low energy Hamiltonian, 

< 'iA,sAVT^,s,>= S _ (23) 

This Rashba like SOC has the same form and sign inde¬ 
pendently on the sublattice where the adatom is placed. 
The Rashba term gets its origin in the broken mirror sym¬ 
metry produced by the adatoms and this is reflected in 
that Ir change sign depending whether the top adatoms 
are located on top or bottom of the graphene layer. 

VI. NUMERICAL RESULTS. 

Adatoms deposited on graphene should be placed at 
minimum energy equilibrium positions. The adsorption 
geometry depends on the particular heavy adatom[22], 
being one of the more interesting that in which the 
adatoms place in hollow positions. When the adatoms 
are intercalated between graphene and the substrate, the 
adatoms form a superlattice commensurated with the 
graphene honeycomb lattice[T8]. It is also plausible to ex¬ 
pect that low energy injected adatoms become deposited 
in random positions. 

In this Section we show numerical results of the elec¬ 
tronic structure of graphene doped with heavy outer 
shell p-orbitals adatoms in three particular cases, i) the 
adatoms are randomly distributed in hollow positions, 
ii) the adatoms form a commensurate supercell with the 
graphene lattice and iii) the adatoms are fully random 
distributed on the graphene sheet. 

In the numerical calculations we consider a periodic 
rectangular graphene supercell of dimensions Lx=NxV^a 


and Ly=Nya, defined by the lattice vectors, A=Nya and 
B = Nx{2h — a). In these expressions Nx and Ny are 
integer numbers. The unit cell contains ANxNy carbon 
atoms located at the graphene lattice positions {Ri+do,}. 
The adatoms are located at positions {r^}. The concen¬ 
tration of adatoms, x, is given by the ratio of number of 
adatoms to the number of carbon atoms. 

The electronic structure is obtained by diagonalizing 
the Hamiltonian, 

H = Ho+ y] {\Z,i,a >< Z,j,a'\+ H.c.) (24) 

where Hq is the pristine graphene Hamiltonian, Eqj^ and 
the second term describes the adatom induced hopping 
between carbon atoms. Because of the periodic bound¬ 
ary conditions, the electronic structure is described using 
the Bloch’s theorem being the electronic states character¬ 
ized by a band index and wave vectors kx and ky that are 
defined in the interval [— p-] and [— p-] respec¬ 
tively. In this geometry and for Ny not being a multiple of 
three, the Dirac cones occur at wavevectors K=(0, 

and K'=(0, Eor Ny multiple of three the two Dirac 

cones overlap at the T point. This overlap does not im¬ 
ply coupling between electronic states in different Dirac 
cones. In order to simplify the analysis of the results, in 
this work we always consider supercells with no overlap¬ 
ping Dirac cones. 

Recent experiments seems to indicate that Pb 
on graphene can induce a large SOC, therefore 
in the numerical calculations we chose Pb as the 
adatom, and we use the tight-binding parame¬ 
ters obtained in reference [18] for Pb atoms on 
graphene, h=0.27nm, Vppa{h)=OAeV, VppT^{h)=-0.6eV, 
ecc=e?/=l-b5eU, ez=1.38eV^ Aso=^AeV and t=2.7eV. 

-- 

•• •• •• •• 
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FIG. 3. (Color online)Graphene supercell with Nx=b and 
Ny=7. Black small circles represent carbon atoms. Larger 
red circles indicate the position of the adatoms. The adatoms 
are located in hollow positions and in this figure we plot a 
particular random realization of disorder. In this figure the 
number of adatoms per carbon atoms is x=17/140. 
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Number of Pb atoms per Carbon atom. 


1 k 

Here ^k=tan“ An adatom located in a hollow po¬ 

sition, at r^, produces a scattering potential of the form, 
s (TzTzdir — Y\). In presence of a density, x, of adatoms in 
hollow positions and neglecting multiple scattering, the 
Green function of the total Hamiltonian is [24] 

G(k, C(j) = Gq( k, C(j) -|- x^sd zGq (k, C(j)t2 ;G(k, cj) . (26) 

Inverting this equation we get, 

f hjjj^x/S.scFz hvpke^^^ \ 
y hvFke~^^^ huo — xAsdz J 


FIG. 4. (Color online) In the inset we plot a band structure of 
a rectangular graphene supercell with x=0.25 Pb atoms per 
carbon. The Pb adatoms induce a gap at the Dirac points 
(0, ^^) and (0, ^^). As commented in the text, the band 
structure is practically independent of the supercell size and 
the gap only depends on the concentration of Pb atoms. In 
the main figure we plot the gap at the Dirac points as function 
of the Pb atoms concentration. The parameters used in the 
calculations are taken from reference m 


A. Adatoms in Hollow positions. 


In this subsection we analyze supercells with differ¬ 
ent sizes and forms and with different concentrations of 
adatoms randomly distributed, but always located in hol¬ 
low positions. In Figj^ we show an example of super¬ 
cell of size Nx=5, with 17 adatoms deposited in 

a random way in hollow positions. In the inset of FigJ^ 
we plot a typical band structure obtained for an adatom 
concentration x=0.25. The adatoms open a gap at the 
Dirac points and, in agreement with the results presented 
in subsection |V A[ the band structure corresponds to a 
Dirac equation in presence of an intrinsic SOC, equations 
Hjandini In our numerical calculations we obtain that for 
atoms adsorbed in hollow positions, the band structure 
always has this form, independently of supercell size and 
form or disorder realization. The SOC depends only on 
the heavy atoms concentration. In Figwe plot the en¬ 
ergy gap as function of the adatom concentration, x. The 
dependence is practically linear and this indicates an al¬ 
most null interference effect between adatoms. 

We understand this linear dependence using Green 
function techniques. The low energy properties of 
an electron are described by the Dirac equation 
hvFcro {skxTx kyTy) and the corresponding Green func¬ 
tion is 


Go(k,u;) 


1 

hw — Hq 
huj‘^ — hv‘^k‘^ 


UJ 

VFke~^^^ 


VFke^^^ 


( 25 ) 


that corresponds to the virtual crystal Hamiltonian 
H=Ho^xAs(JzTz^ describibg graphene in presence of an 
intrinsic SOC of magnitude xA. 




Number of Pb atoms per Carbon atom 


FIG. 5. (Color online) (a) Graphene supercell with 
Lcc=10\/3a and Ly=5a. Black small circles correspond to the 
positions of the carbon atoms. Larger red circles represent 
adatoms. The adatoms form a rectangular lattice commensu¬ 
rate with the graphene supercell. In the figure the dimensions 
of the adatom cell are lx=Lxl9 and ly=Lyl3. (b) Inset, band 
structure for the geometry presented in (a), (b) Main fig¬ 

ure, dependence of the energy gap, as defined in the inset, as 
function of the Pb atoms concentration. Black points are ob¬ 
tained in the graphene supercell fixing the dimension ly—LyjZ 
and changing lx. Red points are obtained fixing lx=Lxl9 and 
changing ly . The parameters used in the calculation are taken 
from reference m- 


to 
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B. Commensurate array of Pb atoms on Graphene. 

Graphene grown on Ir(lll) forms a 9.3x9.3 moire su¬ 
perstructure with a ^25A periodicity [25]. When Pb 
atoms are intercalated under the graphene monolayer, 
the Pb atoms form a rectangular lattice commensurate 
with Ir. Therefore the honeycomb graphene lattice and 
the array of Pb atoms commensurate in a large moire 
super cell [18]. 

In this subsection we study the spin-orbit effects in¬ 
duced by a rectangular array of Pb atoms of dimensions 
lx X ly commensurate with a large rectangular graphene 
supercell of dimensions LxXLy^ see Fig ^a). In the inset 
of Figj^b) we plot, for the geometry shown in FigJ^a), 
the electronic states obtained with the tight-binding pa¬ 
rameters corresponding to Pb. The band structure coin¬ 
cides with the eigenvalues of the Dirac equation in pres¬ 
ence of a Rashba like spin-orbit coupling. The intensity 
of the SO coupling is proportional to the energy gap be¬ 
tween the second conduction band and the second valence 
band. We find that this gap increases linearly with the Pb 
concentration,^, and only depends on the concentration 
of Pb atoms being independent of geometrical details. In 
FigJ^b) we plot the energy gap as function of x for a 
graphene supercell characterized by Nx=10 and Ny=5^ 
and different combination of ix and iy. We obtain the 
same linear dependence for larger graphene supercells. 

These results indicate that Rashba coupling induced 
from different adatoms do not interfere and the total 
Rashba coupling is just the sum of the different contribu¬ 
tions. The Rashba SO coupling induced by Pb adatoms 
have always the same sign, dictated by the broken mir¬ 
ror symmetry, and the linear behavior reveals that in the 
commensurate phase, the adatoms average all possible 
locations in the graphene unit cell. On the contrary, the 
absence of a gap at the Dirac points indicates that the 
contribution to intrinsic SOC from adatoms located in 
different places sums zero. This occurs because the sign 
of the intrinsic SOC induced by adatoms depends on its 
location. A particular example is the case of the intrinsic 
SOC induced by adatoms in hollow position, Eqj^ that 
has opposite sign than the induced by adatoms located 
in top positions Eqj^ 

C. Random Positions 

Einally we compute the SO induced in graphene by a 
concentration of Pb adatoms randomly distributed. We 
study large graphene unit cells, EigJ^a), with different 
concentration of adatoms. The main results of the simu¬ 
lations are that there is not band gap at the Dirac points 
of the band structure, Fig§b), and the Rashba like SOC 
increases linearly with the concentration of Pb atoms, 
Fig|§ These results are independent on the disorder re¬ 
alization and size and form of the unit cell. The depen¬ 
dence of the SOC on Pb concentration is the same than 
in the case of commensurate supercell. This, and the ab- 



Number of Pb atoms per Carbon atom 


FIG. 6. (Color online) (a) Graphene supercell with Lx=7V3a 
and Ly=13a. Black small circles correspond to the positions 
of the carbon atoms. Larger red circles represent adatoms. 
The concentration of adatoms is x=0.2 and they are located 
in a random way. (b) Inset, band structure for the geometry 
presented in (a). Main figure, dependence of the energy gap, 
as defined in the inset, as function of the Pb atoms concen¬ 
tration. The energy gap only depends on Pb concentration 
and is independent on geometry or disorder realization. The 
parameters used in the calculation are taken from reference 

m- 

sence of intrinsic like SOC indicate that in both cases, 
commensurate order and random positions, there is not 
interference effects between adatoms, and the resulting 
SOC is just the sum of the contributions from adatoms 
placed in different positions. 

VII. TOPOLOGICAL PROPERTIES. 

In the previous Sections we have obtained that a 
graphene layer doped with adatoms placed in hollow po¬ 
sitions has a gapped energy band structure similar to that 
obtained from the Dirac equation with an intrinsic SOC. 
On the contrary adatoms randomly or commensurately 
distributed on graphene generate a gapless band struc¬ 
ture that remind that of graphene with Rashba SOC. In 
this Section we check that both adsorption geometries, 
hollow and random, have the same topological proper¬ 
ties that the Dirac equation plus intrinsic and Rashba 
SOC respectively. In order to know the topological prop¬ 
erties, we compute the spin resolved Hall conductivity of 
the system, 

Im n]^\vxPa\'^'^ >< n']<.\vyPcr\n]<. 


(28) 

where |nk > and are the eigenfunction and eigenvec¬ 
tors respectively of the supercell Hamiltonian Eq|^ in 
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tivity is zero. However, in references El and na, it 
was proposed that an exchange field applied to graphene 
in presence of Rashba SOC should open a gap and the 
system would be a non trivial insulator characterized by 
an anomalous quantized Hall effect. We have applied an 
uniform exchange field to the randomly doped graphene 
and we have obtained a gapped band structure, Fig^c) 
and a finite Hall conductivity. 


cr. 


xy — 


'xy 


xy 



adatoms in random positions. 


that proves that adatoms placed randomly on graphene 
generated a Rashba SOC. In this case, and because the 
form of the bands Figj^c), the main contributions to 
the Hall conductivity come from annulus centered at the 
Dirac points, Figj^d). 

The realization of the quantum spin Hall effect and 
quantum anomalous Hall effect for hollow and random 
adatoms respectively , shows the stability of these topo¬ 
logical phases even for a non uniform distribution of the 
adatoms. 


FIG. 7. (Color online)(a)Band structure of a graphene super¬ 
cell (W=5, Ny=7) with a concentration x=0.25 of adatoms 
randomly located in hollow positions, (b) Partial section of 
the Brillouin zone showing in green, the regions that con¬ 
tribute to the spin up Hall conductivity a^y for the case (a), 
(c) Same than (a) for a fully random distribution of adatoms 
and an exchange field of ISmeV. (d) Partial section of the 
Brillouin zone showing in green the regions that contribute to 
the total Hall conductivity a^y + a^y for the case (c). 

the sum the index n and n' run over occupied and empty 
states respectively, is the velocity operator in 

the z/-direction and Per projects the wave function in the 
subspace of spin a. 

In the case of adatoms in hollow positions we have 
obtained 

^xy ~ adatoms in hollow positions. 

that corresponds to a quantum spin Hall system[7]. The 
total Hall conductivity sums zero, as it should be in a 
system with time reversal symmetry. The main contribu¬ 
tions to the Hall conductivity come from circular regions 
centered at the Dirac points, Figj^a)-(b). 

Adatoms placed randomly on graphene do not gener¬ 
ate a gap in the band structure, and the Hall conduc- 


VIII. SUMMARY. 

In this work we have studied the spin orbit coupling 
induced in graphene by heavy adatoms with active elec¬ 
trons residing in p-orbitals. Depending on the location of 
the adatoms we find different induced SOC’s. Adatoms 
located in hollow positions open a gap at the Dirac points 
and induce an intrinsic like SOC. However adatoms ran¬ 
domly placed or commensurate with the graphene lattice 
maintain the system gapless and induce a Rashba like 
SOC. The adatoms only perturb the pristine graphene 
band structure near the Dirac points. 

We find that the SOC induced by the adatoms is ad¬ 
ditive and there is not interference effects or multiple 
scattering. The topological properties of graphene with 
hollow or random adatoms are the same than those of 
the Dirac Hamiltonian in presence of intrinsic or Rashba 
SOC respectively. The finite value of the Hall conduc¬ 
tivity of graphene doped in different geometries indicates 
the robustness of the topological phases against a non 
uniform distribution of the spin orbit coupling. 
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